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Abstract. We develop a complete theoretical description of photoassociative 
Stimulated Raman Adiabatic Passage (STIRAP) near a Feshbach resonance in a 
thermal atomic gas. We show that it is possible to use low intensity laser pulses 
to directly excite the continuum at a Feshbach resonance and transfer nearly the 
entire atomic population to the lowest rovibrational level in the molecular ground 
state. In case of a broad resonance, commonly found in several diatomic alkali 
molecules, our model predicts a transfer efficiency rj up to 97% for a given atom 
pair, and up to 70% when averaged over an atomic ensemble. The laser intensities 
and pulse durations needed for optimal transfer are 10'^ — 10'^ MV/cva^ and several 
fis. Such efficiency compares to or surpasses currently available techniques for 
creating stable diatomic molecules, and the versatility of this approach simplifies 
its potential use for many molecular species. 
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1. Introduction 

The realization of rovibrationally stable dense samples of ultracold diatomic molecules 
remains one of the major goals in the field of atomic and molecular physics. While 
cooling diatomic alkali molecules was seen as a logical next step following the optical 
cooling of atoms, many of the possible applications currently under investigation 
extend beyond atomic and molecular physics. Testing fundamental symmetries based 
on high-precision spectroscopy of ultracold molecules [3 [31 [3] or the attempts to detect 
the time variation of fundamental constants |3] are examples of such applications. 
Another one is ultracold chemistry, where the interacting species and products are in 
a coherent quantum superposition state and could be realized by controlling reactive 
coUisional processes [5]. Important insights about new phases of matter could be 
gained from strong anisotropic dipole-dipole interaction between ultracold dipolar 
molecules [6]. Finally, ultracold polar molecules could also represent an attractive 
platform for quantum computation [71 [8] . Many of those applications require dense 
samples of ultracold polar molecules in the lowest rovibrational state that makes them 
collisionally stable and long-lived. 

Translationally ultracold (100 nK - 1 mK) molecules are produced from an 
ultracold atomic gas by photoassociation (PA) [9j or magnetoassociation (MA) [10]. 
In a typical PA scheme, a pair of colliding atoms is photoassociated into a bound 
electronically excited molecular state that spontaneously decays, forming molecules in 
the electronic ground state. In magnetoassociation, a magnetic field is adiabatically 
swept across a Feshbach resonance, converting two atoms in a matching scattering 
state into a molecule. Both techniques produce weakly bound molecules in highly 
excited vibrational states of the ground electronic potential. Such molecules have to 
be rapidly transferred to deeply bound vibrational states before they are lost from the 
trap due to inelastic collisions. 

Stimulated Raman Adiabatic Passage (STIRAP) [12^ has recently attracted 
significant interest as an efficient way to produce deeply bound molecules, starting 
from Feshbach molecules [131 HI]- It allows to realize high transfer efficiency and 
preserve the high phase-space density of an initial atomic gas. In STIRAP, the laser 
pulses, coupling an initial and a final state to an intermediate excited state, are applied 
in a counter-intuitive sequence where a pump pulse is preceeded by a Stokes pulse. 
During the transfer, the system stays in a "dark" state, i.e., a coherent superposition 
of initial and final states, preventing any losses that would otherwise occur from the 
excited state. By adiabatically changing amplitudes of the laser pulses, the "dark" 
state evolves from the initial to the final state, resulting in nearly 100% transfer 
efficiency [IT. 

Efficient adiabatic passage from the continuum requires laser pulses shorter than 
the coherence time of the continuum [HI [TBI Hi]- The adiabaticity condition of 
STIRAP, r^Ttr ^ 1, where rtr is the transfer time, therefore implies a large effective 
Rabi frequency O for the pulses. In addition, dipole matrix elements between the 
continuum and the bound state are usually small, and so the pump pulse that couples 
the continuum and the excited state would require a very high intensity, which proves 
impractical. Thus the previous STIRAP experiments [I3j . being restricted by the very 
short coherence time of the continuum, used a Feshbach molecular state as an initial 
state. 

The small continuum-bound dipole matrix elements can be dramatically increased 
by photoassociating atoms in the vicinity of a Feshbach resonance. It has been shown. 
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both theoretically and experimentally, that the photoassociation rate increases in the 
presence of a Feshbach resonance by several orders of magnitude [THl UHl 123 1^ - 
This can be explained by considering that delocalized scattering states acquire some 
bound-state character due to admixture of a bound level associated with a closed 
channel, resulting in a large increase of the Franck- Condon factor between the 
initial scattering state and the final excited state. The recently proposed Feshbach 
Optimized Photoassociation (FOPA) technique [21] relies on this enhancement to 
directly reach deeply bound ground state vibrational levels from the scattering 
continuum. Consequently, photoassociation in the vicinity of a Feshbach resonance is 
expected to increase molecular formation rate up to 10® molecules/s [21] . 

In the present work, we combine the approach used in FOPA with STIRAP for 
reducing the required pulse intensity. We predict highly efficient transfer of an entire 
atomic ensemble into the lowest rovibrational level in the molecular ground state. 

The paper is organized as follows. In Section II, we derive a theoretical model 
of a combined atomic and molecular system. Fano theory is used to describe the 
interaction of a bound molecular state with the scattering continuum, represented as 
closed and open channel, respectively. The resulting continuum states are coupled by 
two laser fields to the vibrational target state in the ground state via the intermediate 
excited molecular electronic vibrational state. In Section III, we present the results 
of numerical solutions of the model for several alkali dimers. We find optimal Rabi 
frequencies and profiles of STIRAP pulses for those systems. Finally, we conclude in 
Section IV. 



We consider a three level system as represented in Figured] The ground level labeled 
|1) is the final product state to which a maximun of population must be transfered. 
Typically, this level will be the lowest virational level {v" — 0, J" = 0) of a ground 
molecular potential. This ground level is coupled to an excited bound level |2) of an 
excited molecular potential via a " Stokes" pulse depicted by the blue down-arrow in 
Figure [TJ This level |2) is itself coupled via a pump pulse (red up- arrow) to an initial 
continuum of unbound scattering states jvl/e) of energies e (grey area in Figure [T]). If 
we denote Ci, C2 and C(e) the time dependent amplitudes associated to the final, 
intermediate, and initial states |1), |2), and |5'e), respectively, then the total wave 
function |$) of the system is given by: 



No restriction applies to the definition of the continuum state as it can be 
associated to either a single-channel or a multi-channel scattering state. In this work, 
we consider the multi-channel case in which a bound level \b) associated to a closed 
channel is embedded in the continuum of scattering states |e') of an open channel. 
When the energy of |e') coincides with that of |&), a so-called Feshbach resonance [22] 
occurs. These are common in binary collisions of alkali atoms due to hyperfine mixing 
and the tuning of the Zeeman interaction by an external magnetic field, hence the 
possibility to control interatomic interactions with a magnetic field. Following the 
Fano theory presented in Ref. [23], the scattering state \'i>^) can be expressed as: 



2. Model 




(1) 




(2) 



Efficient formation of ground state ultracold molecules via STIRAP from the continuum at a Feshbach resonanceA 




Figure 1. Schematics: a population from tlie initial state \9e) is transferred to a 
final target state jl) via an intermediate state |2). Both |^e) and |1) are coupled 
to |2) by a pump and a Stokes pulse, respectively labeled Op and f2g. A bound 
level |6) corresponding to a closed channel can be imbedded in the continuum. 



with 



and 



, / /X 1 /r(e') sin A , ^ , ,, , 

6(.,0 = -^^^-cosA^(.-0. (4) 

Here, A = — arct an (2^737^) is the phase shift due to the interaction between \b) and 
the scattering state |e) of the open channel. We assume A S [— 7r/2, 7r/2]. The width 
of the Feshbach resonance, T = 27r|l/(e)p, is weakly dependent on the energy, while 
V^(e) is the interaction strength between the open and closed channels. The position 

of the resonance, ep — Ei, + P J ^^^-73^-^, includes an interaction induced shift from 
the energy of the bound state Ef,. 

If we label Ei the energy of the state \i), the total Hamiltonian H is given by: 

H=J2 E^\i){^\ + /"dee|^',)(^',| + yiight . (5) 

1=1,2 

The light-matter interaction Hamiltonian Viight takes the form: 

Might = -M21 • £s\2){l\ ~ Jde ^2*. • 4 |2> + H.c. , (6) 

where £p g — ep^s^p,S exp(— io^p st) + c.c. are the pump and Stokes laser fields of 
polarization ep,s, respectively, while ^21 and ^2*^ are the dipole transition moments 
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between the states |2) and |1), and |2) and respectively. In this form the 

Hamiltonian aheady takes into account mixing between the bound state of the 
closed channel and scattering states of the open channel. The Schrodinger equation 
describing STIRAP conversion of two atoms into a molecule is: 

ih-^ ^E^C^-fi*^-£*sC2, (7) 
^^~Qf = E2C2- ^21 ■ Es Ci- j de ^2*, ■ Ep C(e), (8) 

ih^^^eC{e)~fll^^-E;C2. (9) 

For simplicity, we set the origin of the energy to be the position of the ground 
state |1), and use the rotating wave approximation with Ci = ci, C2 = 026""^-^*, and 
C(e) — c(e)e^*'^'^^^'^^^*. The Schrodinger equation becomes: 

= -^SC2, (10) 

at 

i—^ = Sc2 - ^sci - de flec{e), (11) 

^^ = A,c{e)~n:c2, (12) 

where S = E2/h— ujs, — e/h— {lus — Wp), and eth is the dissociation energy of the 
ground electronic potential with respect to the state |1). The Rabi frequencies of the 
fields are fls = /t2i • (^s^sl^ (assumed real), Jig = /t2*, • ep£p/h. 

The previous system of three equations can be reduced into a two-equation system 
by eliminating the continuum amplitude c(e) in Eq. (fT2|) . Introducing a solution in the 
form of c(e) = s(e) exp {—iAj) into Eq. (fT2|) . we get 

s = i f dt' n*{t')c2{t')e'^''' +Se{t^O), (13) 

where t = is some moment before the collision of the two atoms. The resulting 
continuum amplitude is 

c^i f dt' r!:(t')c2(i')e'^^^*'"'^ + se(t = 0)e-*^^*. (14) 



Inserting this result into Eq. (Ilip . we obtain a final system of equations for the 
amplitudes of the bound states: 

= - ^sC2, (15) 

ot 

= Sc2 - nsci +1 I de nJt) f dt' nAt')*C2{t')e'^^^''-'^ 
9t Je,^ Jo 

POO 

- / de n,{t)s,{t ^ Q)e~''^^\ (16) 
= 6c2 — ^sci + T — S. 



The third term of Eq. labelled T, corresponds to the back-stimulation 

term, whereas the last term, labelled S, corresponds to the source function. In 
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this source term, the initial ampUtude of the continuum wave function s^{t — 0) 
describing a coUision at to of two atoms with relative energy eo has been discussed in 
various contributions [TH [ini [ITj . A Gaussian wavepacket provides the most classical 
description of a two-atom collision characterized by a minimal uncertainty relation 
between the energy bandwidth of the wavepacket and the duration of the collision: 

(17) 

Futhermore, the Rabi frequency of the field coupling continuum states to 
the state |2) is given by [23] 

iie = ^ —p= sgn(e - ep), (18) 

h ^(r/2)2 + (e-e^)2 

where fl2e is the dipole matrix element between an unperturbed scattering state |e) 
and the state |2), and q is the Fano parameter, expressed as: 



(a*26 ■ep)+Pj 



q= , (19) 

7rF*(e)(^2e • Cp) 

where Cp is the polarization vector of the pump field, and fl2b is the dipole matrix 
element between bound states |2) and \b). The q factor is essentially the ratio of 
the dipole matrix elements from the state |2) to the bound state \b) (modified by the 
continuum) and to an unperturbed continuum state |e). This factor can be made much 
larger than unity, and as will be shown below, the total dipole matrix element from 
the continuum can be enhanced by this factor in the presence of the resonance. The 
magnitude of q can be controlled by the choice of the vibrational state |2). Selecting 
a tightly bound excited vibrational state will increase the bound-bound and decrease 
the continuum-bound dipole matrix elements, resulting in larger q. On the contrary, 
choosing a highly excited state close to a dissociation threshold decreases q. 

Using the expressions given in Eqs. p?|) . (fT8)l . and (fT9|) for the initial amplitude of 
the continuum wave function, the Rabi frequency between the continuum state 
and the excited bound state |2), and the Fano parameter, respectively, we obtain the 
following complete expression for the source term: 

r°° )2 i(,-,o)to 

S^So de g{q,e) sgn{e ~ eF)e ^ e-'^'\ (20) 

with Sq = /i2e • ep£'p/fi,(7r(5e)^/'*, and where the function g{q^ e) is defined as 

5(<?, e) = ; (21) 

We assume that the unperturbed continuum is structureless and the coresponding Rabi 
frequency fl2t ■ BpEp/h depends only weakly on the energy. We also extend tth to — cxd 
to have the initial continuum wavefunction normalized to unity: de |C(e)p = 1. 

We can as well obtain a complete expression for the back-stimulation term T . We 
have: 

2 



T ■ 



fJ-2cep 



£p{t) deg^{q,e) dt' C2it')£pit')e'^^^''-*\ (22) 
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Table 1. Parameters of the Stokes and pump photoassociating pulses providing 
optimal population transfer shown in Fig[2] We use q = 10, 7 = 10* s~^, and 
= ^21 = 0.1 D (1 0=10"^^® esu cm = 0.3934 eag). Rabi frequencies are 
modeled by Gaussians iis,p = ^5 p ^^p (— (t — to ± TS,p))/Tg p, where ± refer to 
the Stokes and pump pulse, respectively. 
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Extending the lower integration limit allows for an analytical solution for the integrals 
over energy and time, leading to the following expression for the back-stimulation 
term: 





2 


h 





TTn£^{t)c2{t) + —iq~l)^£p{t) 



(23) 



3. Results 



In this work, we consider two different cases: first, when T ^ Se, i.e., when the width 
r of the Feshbach resonance is much larger than the thermal energy spread of the 
colliding atoms, and second when T ^ S^. By considering these two limiting cases 
of broad and narrow resonances, more practical expressions for both the source term 
S and the back stimulation term T can be found. The derivation of the final system 
of equations is given in [Appendix A| Here, we only describe the solutions of these 
systems for both broad and narrow resonances. 

Using the parameters of the Stokes and pump photoassociating pulses listed in 
Table [1] for a broad (F = 1 mK) and a narrow (F = 1 /iK) Feshbach resonance, we 
obtain the results depicted in Fig. [21 with the left column corresponding to the broad 
resonance, and the right column to the narrow resonance. The top row shows the 
variation of the Rabi frequencies over the time period required for the population 
transfer calculated using Eqs. (|15p and (|16p along with population in the intermediate 
state 1 2) (middle row) and final state |1) (bottom row). 

For the broad case, we considered a Feshbach resonance with a width F = 1 
mK, which is a typical value for broad resonances (see examples in Appendix A.l), 
and a thermal atomic ensemble with an energy bandwidth = 10 fiK. We see that 
the transfer can reach ^ 97% of the continuum state into the target state |1) (see 
Fig. [2] c). The parameters of the Gaussian laser pulses we used (optimized Rabi 
frequencies, durations and delays of laser pulses) are given in Table [T] the peak 
intensities of the Stokes and pump fields were calculated from Rabi frequencies as 
Ig = cSySn = c{n°h)^/8Trnl^ and Ip = cS^/Sn = c{n°p)^Sj32Tr^/^fil^, where we use 
Eq. (|19p to estimate the continuum-bound dipole matrix element fj,2e ~ t^2b/l'^y{^) = 
V2ii2b/qVnr, resulting in Ip = q^c{n°)^S,T /GAy/^iili,. 
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Figure 2. Time-dependence of tlie Stokes and pump pulses (top row) and 
population in state |2} (middle row) and target state |1) (bottom row) for the 
STIRAP sequence. The left column is for a broad Feshbach resonance, while 
the right column is for a narrow resonance (see Table [l] for values of parameters 
used). The dashed blue lines in the left column are the results obtained without 
resonance, when the parameters are adjusted to obtain the same overall transfer 
fraction as for the broad resonance. Stokes Rabi frequency is in units of 10® s~^, 
while the pump Rabi frequency is in dimensionless units (IGTr/Sg)^''* /t2ee'p£p in 
the broad resonance limit and {2tt /T)^^'^ fi2cep£p in the narrow resonance limit. 
Note that the scale for the Rabi frequencies in the narrow resonance case is 40 
times the scale for the broad resonance, and the magnitude of the pump Rabi 
frequency is enlarged 10 times for better visibility. 



When comparing the results for a broad resonance to the unperturbed continuum 
(i.e., far from the resonance), we find that the source term S is enhanced by the factor 
g{q,ea) (see Eq. (jA.l[) in [Appendix A[ ): 

9(q,<^o) = —f =■ (24) 

This factor has a maximum at 2(eo — ^f)/^ = ^lli with the corresponding maximum 
value vl + g^ « q for q ^ \: hence, the source amplitude is enhanced q times. 
In this limit, all populated continuum states experience the same transition dipole 
matrix element enhancement factor to the state |2), so that the system essentially 
reduces to the case of a flat continuum with an uniformly enhanced transition dipole 
matrix element. One thus expects that in this limit, the adiabatic passage should 
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be efficient, requiring less pump laser intensity when compared to the unperturbed 
{i.e. without resonance) scattering continuum. This is clearly demonstrated in Fig. [5] 
(left column, dashed lines): to reach the same ^ 97% transfer efficiency achieved 
with the broad resonance, a very large pump laser intensity is required if there is 
no resonance in the continuum (Fig. [2] a), while the Stoke laser intensity is basically 
the same. So, the comparison of Rabi frequencies for the broad resonance and no 
resonance cases shows that, to achieve the same transfer efficiencies, the required 
peak pump pulse intensity is about 100 times larger without resonance. Condering 
the intensity used in this particular example, this would lead to intensities in the 
range of 5 x 10^ W/cm^, making STIRAP from the continuum technically impossible 
to achieve without a resonance. This is consistent with the analysis of photoassociative 
adiabatic passage from an unstructured continuum |17j , and the above prediction that 
in the presence of a wide resonance the required pump laser intensity is reduced by a 
factor of '--^ l/?^- 

Results of adiabatic passage in a narrow resonance limit are shown in Fig. [5] 
(right column). We considered a typical value of F = 1 /xK for a narrow resonance 
(see examples in Appendix A. 2) and the ensemble energy bandwidth 5^ = 100 /iK. 
Again, we give the parameters providing the optimal transfer in Table [1] In this 
limit, the transfer efficiency is lower: in the specific case analyzed here, it does not 
exceed 47%. The reason for this lower efficiency is destructive quantum interference 
which leads to electromagnetically induced transparency [251 in the transition from the 
continuum to the excited state. It can be explained using the following argument (see 
Fig. [3]). The limit of a narrow Feshbach resonance corresponds to a weak coupling 
between the bound Feshbach state and the scattering continuum, and thus can be 
neglected in this simplified explanation. The system then can be viewed as consisting 
of bound and continuum states \b) and |c) having the same energy, which are coupled 
by the pump field to a molecular state |2), itself coupled to the state |1) by the Stokes 
field. Assuming that initially all the population is in the state |c), due to the small 
interaction strength between \h) and |c), we can eliminate the state taking into 
account its coupling to |2) by the pump laser as the formation of "dressed" states 
|±) — (|2) ± \h))/\/2. If the dipole matrix element of the \h) |2) transition is 
much larger than that of the |c) |2) transition, the detuning of the "dressed" states 
|A±| = ^ fip'^, Og. As a result, the one-photon couphng of |c) to the excited 
state, as well as two-photon coupling to |1) vanishes, preventing the adiabatic transfer. 
This mechanism is similar to the Fano interference effect, the difference is that the 
continuum is initially populated. One can therefore view it as an inverse Fano effect. 
The effective dipole matrix element of the |c) |2) transition is /i2c ^ l^2b/<lV^- In 
the case we analyzed, q = 10, = F/ \/2(5e = 0.007, and /i2c ~ /^26, which gives ~ 50% 
transfer efficiency. 

The transfer efficiency increases if the Feshbach state is far detuned from the 
populated continuum. Our calculations show that for a Feshbach state detuning 
^ |f^26p/7, the transfer efficiency reaches 70% using the laser pulse parameters in 
Table [TJ We note that the smaller intensity of the pump pulse used for the narrow 
resonance, as compared to the broad resonance, is due to the fact that we used the 
same q = 10 and assumed ^2b = 0.1 D for both resonances. From the definition 
of q, it means that the continuum-bound dipole matrix element /i2e is higher in the 
narrow than in the broad resonance case we considered. This explains the smaller 
resulting pump pulse intensity. The overall conclusion for a narrow resonance is that, 
as opposed to a broad resonance, the presence of the Feshbach resonance prevents 
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|i> |i> 

(a) (b) 



Figure 3. Illustration of the reduction of STIRAP transfer efficiency due to 
destructive quantum interference for a narrow resonance: (a) a simplified level 
scheme where the scattering continuum is modeled by a single state |c) and the 
interaction between the continuum and the Feshbach state \b) is neglected; (b) 
an equivalent scheme, where the strong coupling between the Feshbach state \b) 
and the excited state |c) by the pump field forms "dressed" states |±). We habe 
the Rabi frequency Q2c = ^'^c/^+ + ^-c/^— ~ 0, since C+c = f^-c and 
A+ = -A_. 



one from realizing high transfer efSciencies. It should be noted, however, that the 
destructive quantum interference effect is based on negligible interaction between the 
Feshbach and continuum states during the transfer time, since rtr < ^ T~^. This 
argument shows that already for F > (J^, there is enough interaction to neutralize the 
effect of destructive interference. Therefore, we expect that the broad resonance limit 
can be extended down to F ~ S^, making it applicable to a wide variety of atomic 
species. 

4. Applications to the conversion of an entire atomic ensemble into a 
ground rovibrational molecule gas 

The results of Fig. [5] are for a pair of atoms having a specific mean collision energy 
eo — h{uJs — ujp). Such situation could be realized in very tight traps, e.g., in tight 
optical lattices. For a system with a wider energy distribution, one would like to find 
an ensemble averaged transfer efficiency, and thus one needs to calculate the transfer 
probability P{eo) — |ci|^ for all eg within the thermal spread of energies, and perform 
the averaging as 

where we assume a Maxwell-Boltzmann energy distribution, the pump laser resonant 
with the center of the distribution at (e) = 3/2fcsT, and set the bandwidth of the 
distribution at Se — \J ((Ae)^) = ^3/2A:bT. The results are shown in FiglH In this 
case, while the maximal transfer efficiency in the broad resonance case is ~ 70%, it 
can be achieved with lower laser intensities than in the case of a pair of atoms of Fig. 

EI 

Given the adiabatic photoassociation probability P(e) for two colliding atoms 
with relative energy e, we can calculate the number of atoms photoassociated during 
the time overlap Ttr of the Stokes and pump pulses. During this time, the atom 
with the energy e = /xt'^/2, where \x is the reduced mass, will collide with atoms in 
the volume vr^^wrtr, where n}? is the collision cross-section. The impact parameter 
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Figure 4. Same as Fig. [2] but for the energy averaged transfer. Tlie parameters 
arc listed in Table [2] 



Table 2. Parameters of the Stokes and pump photoassociating pulses providing 
optimal population transfer shown in Fig|4] for averaging over a Maxwell- 



Boltzmann distribution of energies. We use q = 10, 7 
^2[, = ^21 = 0.1 D (1 D=10-^* esu cm = 0.3934 eao). 
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for the collision corresponding to a partial wave with angular momentum J is 
b = {J + l/2)h/p — {J + l/2)h/^/2jIe. The number of collisions that atoms with 
a relative energy in the interval (e, e + de) will experience during the transfer time is 
therefore N{e)de = ■nb'^vTt^p{t)de, where p{e) = 2/9exp (-e/A;BT)^/0r(fcBT)^/^ is 
the spectral density of the atoms {p is the density of the sample) . Finally, J = for 
ultracold s-wave collisions, and the fraction of atoms in the energy interval (e, e + de) 
photoassociated by the two pulses is /(e) = P{e)N{e), or 

- ^r,!f?W2 "trpP(6) exp i-e/ksT). (26) 
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The total fraction of atoms photoassociated by a pair of pulses is / = de /(e) w 
(Pavg)pV27rTtrS^/4/i'^/^-\//cBT, where we assumed that P(e) does not significantly vary 
within the ensemble, and approximated it by the averaged value. Considering as an 
example ^Li atoms at T=100 /iK with an atomic density p — 10^^ cm^'^, an overlap 
time Ttr ~ 1 fj,s, and assuming Pavg — 0.7, the fraction of atoms photoassociated by 
the Stokes and pump pulses is / ~ 2.5 x 10"^: for heavier atoms / ^ 10~^ — 10~^. It 
will therefore require ^ 10* — 10^ pairs of pulses to convert an entire atomic ensemble 
into deeply bound molecules. 

Since only a small fraction of atoms can be transferred to |1) by a pair of 
STIRAP pulses, a train of pulse pairs can be applied to photoassociate the entire 
atomic ensemble. To prevent excitation of molecules in |1) by subsequent pulses, 
they have to be removed before the next pair of pulses is applied. This could be 
realized by applying, after each pair of Stokes and pump pulses, a relatively long 
pulse resonant to a transition from |1) to some other vibrational level in the excited 
electronic potential which decays spontaneously to a deep vibrational state in the 
ground electronic potential. This long pulse would optically pump molecules out of the 
state |1) to deeper vibrational states in the ground electronic potential. It therefore has 
to be longer than the spontaneous decay time of the excited state. Care has to be taken 
that the excited state does not decay back into the scattering continuum. This would 
empty the |1) state and deposit molecules into ground potential vibrational states 
according to Franck-Condon factors before the next pair of pulses arrives. Finally, 
after all atoms have been converted into molecules the recently demonstrated optical 
pumping for molecules method [26] can be applied, which would transfer molecules 
from all populated vibrational states into the ground level v = 0. 

The optimal strategy is to actually choose an excited state that decays mostly 
to the V = level. This would allow one to avoid storing molecules in unstable 
vibrational states and using the optical pumping method. If such a state cannot be 
directly reached from a four-photon STIRAP transfer can be applied [27], which 
provides efficient transfer to deeply bound molecular states. It allows one to choose 
the final state from which the excited state decaying predominantly to w = can 
be easily reached. In this case rotational selectivity can also be preserved, since only 
V = 0, J = and v = 0, J — 2 states will be populated. 

The total time required to photoassociate the whole atomic ensemble and transfer 
it to the V = level can be estimated as follows. As the numerical results show, 
adiabatic passage requires ~ 5 fj,s, the follow-up pulse emptying state |1) can have a 
^ 100 ns duration, if the excited state lifetime is tens of ns, resulting in the whole 
sequence ~ 6 fxs. Then the train of 10* — 10^ pulse pairs will take ~ 0.1 — 10 s. The 
final step, optical pumping to the v — level, requires ~ hundred /iS, so the overall 
formation time is ~ 0.1 — 10 s. Given an illuminated volume ~ 1 mm'^ and an atomic 
density p ^ 10^^ cm^'^ the resulting production rate is expected to be 10^ — 10^" 
molecules/s. This compares well with the recent experiment on STIRAP production 
of ground state KRb molecules starting from the Feshbach state, where the entire 
cycle including creation of Feshbach molecules takes ~ 10 s [13] . 

5. Conclusion 

Combining photoassociation and coherent optical transfer to molecular ground 
vibrational states can allow one to convert an entire atomic ensemble into deeply 
bound molecules, and to produce a high phase-space density ultracold molecular gas. 
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We have analyzed photoassociative adiabatic passage in a thermal ultracold atomic 
gas near a Feshbach resonance. The presence of a bound state imbedded in and 
resonant with scattering continuum states strongly enhances the continuum-bound 
transition dipole matrix element to an excited electronic state, thus requiring less 
laser intensity for efficient transfer. In the limit of a wide resonance when compared 
to the thermal spread of collision energies, the dipole matrix element is enhanced by 
the Fano parameter q. Choosing a tightly bound excited vibrational state, q can be 
made much larger than unity, resulting in the intensity of the pump pulse required 
for efficient adiabatic passage to be ~ times smaller than in the absence of the 
resonance. We modeled the adiabatic passage using typical parameters of alkali dimers 
and found intensities and durations of STIRAP pulses providing optimal transfer. 
Intensities of the pump pulse, coupling the continuum to an excited state, were found 
to be a few kW/cm^, which is ~ 100 times smaller than without resonance. Optimal 
pulse durations are several /is, resulting in energies per pulse ^ 10 /iJ for a focus area 
of 1 mm^. 

If the Feshbach resonance is narrow compared to the thermal energy spread of 
colliding atoms, adiabatic passage is hindered by destructive quantum interference. 
The reason is that electromagnetically induced transparency significantly reduces 
the transition dipole matrix element from the scattering continuum to an excited 
state in the presence of the bound Feshbach state. In the narrow resonance limit, 
photoassociative adiabatic passage is therefore more efficient if the resonance is far- 
detuned. 

Due to low atomic collision rates at ultracold temperatures, only a small fraction 
of atoms can be converted into molecules by a pair of photoassociative pulses. To 
convert an entire atomic ensemble, a train of pulse pairs can be applied. We estimate 
that lO"' — 10^ pulse pairs will associate an atomic gas of alkali dimers with a density 
10^^ cm~^ in an illuminated volume of 1 mm'^ in 0.1 — 10 s, resulting in extremely high 
production rates of 10* — 10^° molecules/s. High transfer efficiencies combined with 
low intensities of adiabatic photoassociative pulses also make the broad resonance limit 
attractive for quantum computation. For example, a scheme proposed in [28] can be 
realized, where qubit states are encoded into a scattering and a bound molecular states 
of polar molecules. To perform one and two-qubit operations, this scheme requires a 
high degree of control over the system, which our model readily offers. 

Finally, marrying FOPA and STIRAP is a very promising avenue to produce large 
amounts of molecules, for a variety of molecular species. In fact, although we described 
here examples based on magnetically induced Feshbach resonances, such resonances 
are extremely common, and can be induced by several interactions, such as external 
electric fields or optical fields. Even in the absence of hyperfine interactions, other 
interactions can provide the necessary coupling, such as in the case of the magnetic 
dipole-dipole interaction in ^^Cr [291130]. 
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Appendix A. Adiabatic passage in the limits of broad and narrow 
Feshbach resonances 

In this appendix, we discuss Eqs. (|15p and (|16p for various relative widtlis of the 
Feshbach resonance F with respect to the thermal energy spread of the colliding 
atoms. We first describe the case of a broad resonance, i.e., when the width of the 
Feshbach resonance greatly exceeds the thermal energy spread (F :§> S^), and second 
consider the opposite situation of a narrow resonance (F <^ d^). Finally, we briefly 
present the case where there is no resonance. 



Appendix A.l. Limit of a broad Feshbach resonance T ^ 5e 



The typical thermal energy spread for colliding atoms in photoassociation experiments 
with non-degenerate gases is <5e ~ 10 — 100 /xK. The broad resonance case occurs 
for resonances having a width of several Gauss 1 niK), for which we have 
T /5e ~ 10 — 100. A wide variety of systems exhibit broad resonances. For instance, 
they can be found in collision of ^Li atoms at 834 G for the \f = 1/2, m/ = 1/2) 
channel (F = 302 G= 40 mK) and in '^Li at 736 G for the \ f = 1,to/ = 1) channel 
(F = 145 G = 19 mK). We note here that these values of F are slightly different than 
the "magnetic" width AB usually given and based on the modelling of the scattering 
length. 

The source function can be readily calculated from Eq.([20ll by noticing that the 
Rabi frequency term can be set at e = eg corresponding to the maximum of the 
Gaussian function in the integrand. Using the function g{q, e) defined in Eq. (|21|) . the 
result takes the form 

Sno-ics9{q, eo)sgn(eo - ep), (A.l) 

where 5'no-rcs is the source function without a resonance given below in Eq. (|A.10|) . 
Strictly speaking, this expression is valid for |e_F — eqI > ^e, but since F ^ (5^ Eq. (jA.l|) 
is a good approximation for a wide range of detunings ep — eo- 

The back-stimulation term (|23p can be further simplified in the limit of a broad 
resonance. In this case, both C2{t) and £p{t) change on a time scale ~ l/<5e, i.e., slowly 
compared to the decay time ~ h/T of the exponent. Therefore, we can rewrite 
as: 





2 






h 





1 



1 + 2i{eF - Hlus - u;p))/r 



C2{t)£l(t) 



The system ((15 



mi 

.dci 
■ 9^ 
.dc2 
' 9^ 



in the case of a broad resonance becomes 

= - VtsC2, 

= - ^sci - Sw + {& - i"f)c2 

- i7rft|il„o-rcs(t)P 



1 + 



1 + 2?:(eF - hiios - LUp))/r 



(A.2) 



(A.3) 



C2,(A.4) 



where ^no-res = {l2eSp£p/h IS the continuum-bound Rabi frequency in the absence 
of resonance. We also added a spontaneous decay term 7C2, assuming that the 
excited molecules dissociate into high energy continuum states and the resulting atoms 
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leave a trap. From Ea. (|A.ip . one can see that in a broad resonance case, the source 
amplitude is enhanced by the factor g{q, eg) = (9 + 2(60 — e_F)/r)/-\/l + 4(eo ~ epY /^'^ 
when compared to the unperturbed continuum case. This factor has a maximum at 
2(eo — e_F)/r — with the corresponding maximum value ^max ~ 9 for g S> 1. 

Appendix A. 2. Limit of a narrow Feshbach resonance F <C 5e 

This situation occurs when the width of the resonance is of the order of a few micro- 
Gauss or less. Examples of narrow resonances include ^Li^'^Na at 746 G for the 
1/1 = 1/2, m/i = 1/2)1/2 = l,mf2 = 1) channel (F = 7.8 niG = 1 fiK) [IJ, or 
^Li^^Rb at 882 G for the |/i 1/2, m/i = l/2)|/2 = l,m/2 = 1) channel (p-wave, 
F = 10 mG = 1.3 ^K). 

We note that the source term expressed in Eq. ([20| can be rewritten in a time 
representation: 



/OO 
e-(^'-^Df{I,{^\r-r,-T'\) 
-OC 



- L_i(C|t - to - r'l) - tqimir - Tq - t'|) 

-Lo^T -TO- r'|))sgn(r -To- r'))dr'] , (A.5) 

where we introduced the dimensionless variables t — tS^/V^h, D — {ep — eo)/V^St, 
^ = r/-\/2<5e; /o,i and io,-i are modified Bessel and Struve functions. One can see 
from this expression that the source function is a sum of the pure source function of 
the unperturbed continuum, given by the first term in square brackets, and of the 
admixed bound state, given by the integral. The coefficient ^ = F/\/25e, which is 
the ratio of the Feshbach resonance width to the width of the thermal energy spread, 
gives the ratio of contributions from the bound state and the unperturbed continuum, 
respectively. 

It is then easier to notice that in the limit of a narrow resonance, the Gaussian 
function in the integrand of Eq. ljA.SP is much narrower than the Bessel and Struve 
functions, which change on the time scale ~ 1/^. Therefore the source term can be 
aproximated as: 

Sn = SoV2^See-'^'">/''-^'^''-'^'''>^'[e-^^-^°'>' 

+ ^V^e2'^-^'(/i(e|r - rol) - L^i{^\t - to\) 

- iqilo^T - rol) - Lo(C|t - ro|))sgn(r - To))]. (A.6) 

Since ^ ^ 1, the real part of the source function is given by the first term in 
the square brackets, which is a pure continuum source function, while the imaginary 
part is due to the admixed bound state and its magnitude depends on the product ^q. 
Using asymptotic expansions of modified Bessel and Struve functions Io{x) — Lo{x) 
—2/ttx, Ii{x) — L-i{x) — *■ —2/ttx'^, it is seen from Ea. (jA.6|) that the contribution to 
the source function from the bound state decays on the time scale |t — ro| ^ 1/^, 
while the contribution from the unperturbed continuum decays on the time scale 
|r-To|~l«l/e 

In the limit of a narrow resonance the system (fT5 |) -(fT6 l) becomes: 

= - f^5C2, (A.7) 
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.dC2 



dt 



VtsCi - Sn + (6 ~ ij)c2 

^ ^ 2 



^h£lc2 + —{q-if£p(t) 



(A.; 



Appendix A. 3. Continuum without resonance 



Finally, let us consider the case of a continuum without resonance. In this case the 
continuum-bound Rabi frequency Eq. (jl8p is: 

(A.9) 
(A.IO) 

(A.11) 



— ^no — ICS — MSe ' ^p/^i 

and the source function is 

The back-stimulation term ([^5]) reduces to 

2 2 
fL2e-ep/h TTh£pC2 = nh\Qno-i-cs{t)\ C2, 



and the system p5|) - (ri6|) takes the simple form: 

dci 

.dc2 



^SC2: 



(A.12) 



dt 



^SCl + {5 ~ i-f)c2~ iTTh\U.r,o-xcs{t)\ C2 - 5'no-rcs- (A. 13) 
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